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Abstract 

A many body Hamiltonian comprising pairing, quadrupole-quadrupole and 

spin-spin interaction is treated within a projected spherical basis with the 

aim of describing the detailed structure of the magnetic states of orbital and 

spin-flip nature in ^^^'^^^~^^^Sm. The mean field for the single particle motion 

takes care of the volume conservation of a phenomenological core which results 

in having a complex dependence on deformation for single particle energies. 

Several experimental features like the total orbital and total spin strengths, 

the shape of the orbital strength distribution, the hump structure of the spin 

strength distribution, the quadratic dependence of the orbital strength on 

nuclear deformation, the saturation effect of the orbital strength with respect 

to the Casten's factor P are very well described within a quasiparticle random 

phase approximation formalism. Comparison of present results with those 

obtained in some previous treatments are also given. 
PACS number(s): 21.10.Re, 21.60.Ev, 27.10.+j, 27.70.+q 
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I. INTRODUCTION 



The collective magnetic dipole mode (Ml) was one of the central subjects in the last 
two decades. The interest has grown with time due to the wealth of experimental data pro- 
vided by several groups performing either (e, e') experiments at backward angles with high 
resolution or (7, 7') measurements using brehmstralung radiation. The theoretical work in 
this field preceded the experiment. Indeed Lo ludice and Palumbo [|T| predicted the ex- 
istence of this state from an anti-phase angular oscillation of proton and neutron axially 
symmetric and rigid systems. The model is known under the name of two rotor model 
(TRM), and the predicted mode as scissors mode. Later on lachello predicted the existence 
of such a mode at lower energy within the Interacting Boson Approximation (IB A) 0. 
The experimental discovery of this state in ^^^Gd, by Achim Richter and his collaborators 
0], in an (e, e') experiment, was a stimulus for both experimental [^-[T0[| and theoretical 
rT|^6[ intensive activities. Thus, many descriptions have been proposed and some authors 



reactivated some old own papers devoted in the past to closely related topics. Both phe- 
nomenological [llI HT3|JT5|JT6|j3^ and microscopic |lT^ , p!7| -p^ descriptions have been proposed. 
The phenomenological formalisms provide an intuitive basis for the description of the new 
mode, although many details such as the fragmentation and the magnitude of the total 
strength cannot be accounted for. That is the reason why several improvements have been 
proposed along the time. The rigid motion of the proton-neutron system was replaced by 



a share motion pO[ or the motion of two liquid drops ||ll| , p!7|| . The Interacting Boson Ap- 
proximation was extended by distinguishing between the proton and neutron bosons ||3^ 



and pointed out that the F spin symmetry plays an important role for the structure of the 
magnetic state. Also the Coherent State Model [^] was generalized to GCSM [jl2|, which 
provides a new boson description of Ml states. The microscopic formalisms use the QRPA 
approach within a deformed single particle basis and predict indeed that the Ml strength is 
distributed among many QRPA states. The nice feature brought by the microscopic theories 
is that the low lying states with relatively large Ml transition probability are excited mostly 



2 



by the orbital part of the transition operator. Such type of states are conventionally called 
scissors like states since their structure confirm the hypothesis of the TRM model that the 
mode is due to the convection current created by the protons motion. 

Early microscopic description of the state as a doorway state orthogonal to the spurious 
rotational modes predicted fragmentation The fragmentation effect was con- 



firmed experimentally and as a consequence a new question was opened, namely, whether 
there are orbital states in the region beyond 4 MeV where the spin-flip excitations dominate. 
A long standing debate was about the mixture in the predicted mode, of spurious components 
of isoscalar nature |^2|. Several projection procedures were proposed [p!^ , ^ , ^ , ^ , ^ ,^,^ 



to purify the theoretical state from spurious components. 

Many review articles discussing the results on this topic appeared in journals or pro- 
ceedings of conferences and summer schools. Since we focus our attention on a specific 
description, we might unavoidably omit some important contributions to the field from our 
list of references and therefore, we advise the reader to consult the most recent reviews from 

Refs. mwn. 



The present paper is devoted to the study of the orbital and spin excitations in the even- 
even isotopes of Sm, using a projected single particle basis [Q. Some ingredients like the 
projected states and the model Hamiltonian are similar to those used in Ref. |^6|. The new 
point of the present investigation consists of the definition of the mean field. Here the mean 
field of the single particle motion is defined by averaging a particle core Hamiltonian on a 
coherent state describing a phenomenological deformed core. The coupling term comprises a 
monopole-monopole and a quadrupole-quadrupole interaction. The monopole term is calcu- 
lated by exploiting the restriction for the volume conservation. In this way the single particle 
energies, approximated by the average of the mean field Hamiltonian on the projected states, 
have a complex dependence on the deformation parameter. We recall that in Ref. the 
single particle energies depend linearly on deformation. Therefore our aim is to investigate 
whether by considering more realistic single particle energies one obtains an improvement 
of the quantitative description given in Ref. fS^. As we shall see at the end of this paper 



many details shown by the experimental data are nicely described by our approach which in 
fact confirm our expectation that in order to describe the non-global features of the system 
one should define properly the mean field of the single particle motion. 

The results are described according to the following plan. In Section 2, we present briefly 
the projected single particle basis. Also, here we describe the mean field and derive analytical 
expressions for the single particle energies. The model Hamiltonian and the QRPA states are 
introduced in Section 3. The reduced probability for exciting these states from the ground 
state is also given. Applications to the even-even isotopes of Sm are described in Section 4. 
The main results are summarized in Section 5 



II. PROJECTED SINGLE PARTICLE BASIS 

Technical details about the construction of the angular momentum projected single par- 
ticle basis appropriate for the description of the single particle motion in a deformed mean 



field generated by the particle core interaction, are given in Ref. [Q. Aiming at a self- 
consistent presentation, we give here the main ingredients of the method used in the above 
quoted reference. 

The single particle mean field is determined by a particle-core Hamiltonian: 

H = Hs„, + H,ore-Mujy J2 E "v^Am- (2.1) 

A=0,2 -A<^t<A 

Here denotes the spherical shell model Hamiltonian while Hcore is a harmonic 

quadrupole boson (6+) Hamiltonian associated to a phenomenological core. The interac- 
tion of the two subsystems is accounted for by the third term of the above equation, written 
in terms of the shape coordinates aoo, a2^f The quadrupole shape coordinates are related to 
the quadrupole boson operators by the canonical transformation: 

«2, = ^(4, + (-rv.), (2.2) 

where k is an arbitrary C number. The monopole shape coordinate is to be determined from 
the volume conservation condition. In the quantized form, the result is: 
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(2.3) 



Averaging H on the eigenstates of Hs^, hereafter denoted by \nljm), one obtains a deformed 
boson Hamiltonian whose ground state is described by a coherent state 

^, = expK6+ -62o)]|0)fe, (2.4) 

with |0);, standing for the vacuum state of the boson operators and d a real parameter which 
simulates the nuclear deformation. On the other hand, the average of H on ^ g is similar 



to the Nilsson Hamiltonian |^9[|. Due to these properties, it is expected that the best trial 



functions to be used to generate, through projection, a spherical basis are: 

= \nljm)mg. (2.5) 

The upper index appearing in the l.h. side of the above equation suggests that the product 
function is associated to the particle core system The projected states are obtained in the 
usual manner by acting on these deformed states with the projection operator 

PiiK = / Di,^*{Q)R{Q)dn. (2.6) 

A certain subset of projected states is orthogonal: 

^if;{d) = Ar^,^pU\nijm,]. (2.7) 

The main properties of these projected spherical states are: a) They are orthogonal with 
respect to I and M quantum numbers, b) Although the projected states are associated 
to the particle-core system, they can be used as a single particle basis. Indeed when a 
matrix element of a particle like operator is calculated, the integration on the core collective 
coordinates is performed first, which results in obtaining a final factorized expression: one 
factor carries the dependence on deformation and one is a spherical shell model matrix 
element, c) The connection between the nuclear deformation and the parameter d entering 
the definition of the coherent state (2.4) is readily obtained by requiring that the strength 



of the particle-core quadrupole-quadrupole interaction be identical to the Nilsson deformed 
term of the mean field: 

Here Q± and denote the frequencies of Nilsson's mean field related to the deformation 



45/167r/3 by: 



= (^)^/^ = (^)-^/^- (2.9) 



The average of the particle-core Hamiltonian H' = H ~ Hcore on the projected spherical 
states defined by Eq.(2.7) has the expression 



90~ 



(2.10) 



Here we used the standard notation for the Clebsch Gordan coefficients C^^j'^^m- ^j'^ stands 
for the following integral 

Ap = Pj{x)\P2{x)fexp\d^P2{x)]dx, A; = 0, 1 (2.11) 

JO 

It is worth mentioning that the norms for the core's projected states as well as the matrix 
elements of any boson operator on these projected states can be fully determined once the 
overlap integrals defined in (2.11), are known [0. Since the core contribution does not 



depend on the quantum numbers of the single particle energy level, it produces a shift for 
all energies and therefore is omitted in Eq.(2.10). However when the ground state energy 
variation against deformation is studied, this term should be considered as well. 

The first term from (2.10) is, of course, the single particle energy for the spherical shell 
model state \nljm). According to our remark concerning the use of the projected spherical 
states for describing the single particle motion, the average values e^;^ may be viewed as 
approximate expressions for the single particle energies in deformed Nilsson orbits. We may 
account for the deviations from the exact eigenvalues by considering, later on, the exact 
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matrix elements of the two body interaction when a specific treatment of the many body 
system is apphed. It is worth mentioning that the dependence of the new single particle 
energies on deformation is similar to that shown by the Nilsson model. This is clearly seen 
in Fig. 1 and Fig. 2 described in section 4. We remark that the energies shown in the above 
mentioned plot depend on deformation in a different manner than those obtained by one of 



the present authors in Ref. Indeed therein, they depend linearly on deformation while 
here non-linear effects are present. The difference between the two sets of energies is caused 
by the fact that here the volume conservation condition was used for the monopole shape 
coordinate while in Ref. ||3^ this term was ignored and moreover there is no account of the 
volume conservation condition. In this respect, one expects that this difference in the single 



particle energies will induce some modifications of the results of Ref. ||3^ concerning the 
fragmentation process of the Ml strength. Although the energy levels are similar to those of 
the Nilsson model, the quantum numbers in the two schemes are different. Indeed here we 
generate from each j a multiplet of (2j + 1) states as I, which plays the role of the Nilsson 
quantum number Q, runs from 1/2 to j. On the other hand, for a given I there are 2/ + 1 
degenerate sub-states while the Nilsson states are only double degenerate. As explained in 
Ref. |PB[, the redundancy problem can be solved by changing the normalization of the model 



functions: 

Q/Af|*&a/Af) — 2. (2-12) 

M 

Due to this normalization, the states ^aiM used to calculate the matrix elements of a given 



operator should be multiplied with the weighting factor y2/(2/ -|- 1). 

Concluding, the projected single particle basis is defined by Eq.(2.7). The projected 
states are thought of as eigenstates of an effective Hamiltonian, with the corresponding 
energies given by Eq.(2.10). 
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III. THE MODEL HAMILTONIAN AND THE Ml STATES 



In the states introduced in the previous Section, we consider a system of Z protons 
and N neutrons interacting among themselves in the following manner. The nucleons of 
similar isotopic charge interact through pairing interaction, dipole-dipole and quadrupole 
quadrupole interaction. Protons and neutrons interact only through dipole and quadrupole 
two body interaction. The many body Hamiltonian describing such a system is written in 
second quantization as: 

aIM;T r ^ 



- ^ E xi'Xir)a^ir') - 1 ^ xi'J,Q,,ir)Qt,iT'). (3.1) 



flTT' fl 



TT' 



Note that although the model Hamiltonian is formally identical to that used in Ref. ^6 
the mean field parts of the two schemes are different from each other. Therefore even if 
the procedures used are similar, the final results are different. Before applying the QRPA 
treatment, it is convenient to rewrite the QQ interaction by re-coupling the initial and final 
nucleons in a different way and keep in the final result only the dipole-dipole term. By doing 
this, the new expression is no longer separable. The many body Hamiltonian is written first 
in quasiparticle representation, by using the Bogoliubov-Valatin transformation, defined by 
the coefficients Uai and Vai- At this stage one can define the QRPA dipole phonon operator 
which creates, by acting on the vacuum state |0), the QRPA states: 

|1V,A) = E [X^,{r)A+,iab-r)-Y^,ir)Ai,iab-T)]\0) ^ r+ (A)|0), (3.2) 

a<b;T 

where y4+(a6; r) and A^{ab; r) denote the dipole operator constructed out of two quasipar- 
ticle creation operators and and two annihilation operators Oq, and a^, respectively 

= E^4"XX«^' M<^b;r) = {At,{ab;r)y {-y-^. (3.3) 

q:/3 

The phonon amplitudes are obtained by solving the QRPA equations with the normalization 
condition: 
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E{i^a.i' = 1- (3.4) 

a<b;T 

The solutions of the QRPA equations are labeled by the index A = 1,2, ... defined by the 
ordering relation for the corresponding energies {ujx): uji < 002 < 003 < ... The QRPA states 
can be excited from the ground state, the vacuum state of the phonon operator, by the 
transition operator: 

= J]^j:[9t"i^)l,i^^'r) + gf^iT)s,it,T)]. (3.5) 

where the effective gyromagnetic factors are obtained by accounting for the renormalization 
induced by the angular momentum carried by the core system: 

gf^ = 9i- 9c, gf^ = 9s- Qc- (3.6) 

The reduced transition probability is defined byQ: 

S(M1;0+ ^ 1+) = B{M1 t) = |(0||A^(M1)||1+,A)|2. (3.7) 

The involved matrix element has the expression: 

(0||A1||1+, A) = \I^Y: i(^7p(r)U(r) + lg:"iT)a^,{T)) 

X (t/a(r)H(r) - U,{T)Va{r))[Xl,{T) - Y^,{r)l I = V2TT1. (3.8) 

with obvious notations for the reduced matrix elements for the orbital angular momentum 
(1) and spin operator (|). 

IV. APPLICATION TO THE EVEN-EVEN SM ISOTOPES 

The QRPA equations were solved for I44,i48-i54gj^ isotopes. The single particle energies 
are those given by Eq. (2.10), with the spherical shell model energies enij determined by the 
frequency ujq and the strengths C and D for the spin orbit and terms, respectively 



^Throughout this paper the Rose convention for the reduced matrix element is used 
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AlA-^/^MeV, C = -2uJoX, D = -uoXf^, (4.1) 



where x = 0.0637, /i = 0.6 for protons and /i = 0.42 for neutrons Il49|j50[| . The single particle 
energies can be considered either as functions of the nuclear deformation 6 by means of the 
frequencies fl±, Q^, or as functions of the deformation parameter d due to Eq. (2.8). In the 
later situation, an additional parameter is introduced, fc, which remained undetermined when 
the shape coordinates were transformed in boson variables. In order to compare the present 
results with those reported in Ref. 0] it is useful to define the strength of the particle core 
qQ interaction so that the corresponding matrix elements have similar expressions in the 
two descriptions. It results 

X.= ^./A. (4.2) 
8fc V 27r ^ ^ 

For illustration, proton and neutron single particle energies are presented as functions of 
the deformation parameter d in the figures 1 and 2, respectively, for ^^^Sm. The particle- 
core coupling constant X^. was taken equal to that given in Ref. |^. In principle the 



deformation parameter should be determined by requiring that the system of Z protons 
and N neutrons filling the energy levels in Figs. 1 and 2 respectively, in accordance to the 
Pauli exclusion principle, has a minimum energy with respect to changing the deformation. 
By this procedure, one obtains the equilibrium value of d specified by the dashed vertical 
line. However in the generalized coherent state model (GCSM) |12|, which is also based on 
coherent states description, the deformation parameter was fixed by a fit of the experimental 
E2 strength to the lowest 2+ states. It turns out that this procedure works well for deformed 
nuclei but not unambiguously for near spherical nuclei. The value obtained for the case of 
^^^Sm, by using this method, is mentioned by vertical full lines in Figs. 1, 2. The deviation 
of this value from the equilibrium deformation might be corrected in two ways; either by 
modifying the strength (or, equivalently, the parameter k) so that the two values are the 
same, or by considering the normalization of the single particle energies due to the motion 
of the core. Indeed, we recall that the term of the model Hamiltonian which corresponds 
to the core motion was set to zero when the average on the coherent state describing the 
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phenomenological core was performed. The average of the core Hamiltonian is a function of 
d which influences, to a large extent, the equations for energy minimum. The Fermi energies 
for protons and neutrons are pointed by arrows in the corresponding graphs. 

The strength parameters of the model Hamiltonian which were used in our calculations 
are collected in table 1. As we have already mentioned, we want to outline the effect on 
the Ml strength coming from the fact that we use a different set of single particle energies 
than in Ref. Such an effect could be best determined if the strength parameters used 
in the two descriptions were close to each other. Indeed the parameters from Table 1 
are, with few exceptions, practically the same as in Ref. [^. We modified the pairing 
strengths so that the gap parameters produced with the new set of single particle energies 
were the same as in the previous treatment. Moreover these are the same as those used by 
Hamamoto and Magnusson in Ref Only for ^^^Sm and ^^^Sm one notices a deviation 
of the proton gap of about lOOkeV. Also the strength of the QQ interaction for ^^''Sm and 
^^^Sm is smaller in the present case by a factor of about 2. The strengths of the pairing 
interaction are consistent with the results obtained by Nilsson while the strengths of 
dipole and quadrupole interactions are close to the parameters used by de Coster and Heyde 
PP|^T| and by Hamamoto and Magnusson . The results of the BCS treatment are shown 
in terms of quasiparticle energies for ^^"^Sm, in Figs. 3 and 4. From these figures, one sees 
that despite the fact that the considered nucleus is a deformed one, the quasiparticle energies 
exhibit a shell like structure manifested by some energy gaps. The quasiparticle energies 
were used to solve the QRPA equations. Further, the transition operator given by Eq. (3.5) 
is used to excite the resulting states. The reduced transition probability for such excitation 
to take place is calculated by means of relation (3.8) which contains in a factorized form the 
contribution coming from the BV and QRPA transformation amplitudes as well as from the 
matrix elements of the orbital and spin operators. 



11 



(MeV) (MeV) (MeV) 
^^^Sm 0.400 2.330 1.500 20.70/A 1.592 19.51/A 0.951 0.797 0.0171 1.196 
i^^Sm 0.500 2.750 2.800 19.66/A 1.328 17.35/A 1.060 0.600 0.0166 1.164 
^5°Sm 1.482 3.770 1.100 20.52/A 1.227 18.16/A 1.331 0.600 0.0164 0.533 
^52Sm 1.850 4.980 0.688 18.27/A 0.865 17.08/A 1.049 0.500 0.0162 0.658 
is^Sm 2.290 5.580 0.611 18.20/A 0.777 15.80/A 0.721 0.850 0.0160 1.118 

Table l.The parameters entering the many body Hamiltonian (3.1) are listed. Gap energies for 
proton and neutron systems are also given. The strength of the QQ interaction is multiplied by 
the fourth power of the oscillator length b{=^/K/JmLJo). 

The results for Ml strength are presented in several figures for the energy regions of orbital 
and spin transitions. The total strengths characterizing the intervals 2-4 MeV and 4-10 MeV are 



listed in Table 2. Comparing the data given there and the corresponding ones given in Ref. |36| one 
immediately concludes that the single particle energies have a notable effect on both the orbital 
and the spin excitations in the even-even Sm isotopes. For the gyromagnetic factors of orbital 
angular momenta, we took their bare values, i.e. gi{p) = 1 and gi{n) = 0. The spin gyromagnetic 
factors are quenched by a factor 0.85, as in Ref. |2^, in order to account for effects determined 
by nucleon excited states which are not described in the present paper. For the gyromagnetic 
factor induced by the core angular momentum, we took the value gc = 0.166 which corresponds 
to a ratio of the renormalized neutron to proton orbital gyromagnetic factors equal to 0.2, as 



suggested by the phenomenological models IBA2 [34| and GCSM [52| . Our analysis suggests 
a small contribution due to the core gyromagnetic factor. For orbital transition strength this 
confirms the scissor character of the transition, i.e. it is described by an operator proportional to 
i9p ~ 9n){ifi — l^)- As far as the spin-flip transitions are concerned, the correction due to gc is small 
since this is small comparing it to the spin gyromagnetic factor. The theoretical and experimental 
spin-spin strength distributions are compared with each other in Fig. 5. It is interesting to remark 
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that all peaks in the two panels showing experimental and theoretical results can be easily put in 
one to one correspondence. The difference between the two pictures consists of that two of our 
predicted peaks at about 7.2 and 7.8 MeV are higher than the corresponding experimental peaks 
7.5 and 7.9MeV. Also in the low part, the experience shows a broad peak of 0.05 ^% while in the 
present calculations one finds a narrow peak which is 0.4 /x^. 





i44Sm 


i^^Sm 


isosm 


i52Sm 


is^Sm 


Si 


0.260 


0.401 


0.954 


2.324 


2.508 


S2 


5.656 


5.473 


7.484 


7.755 


7.961 



Table 2. Predictions for the total orbital strength, below 4 MeV, and total spin-flip transitions 
to the states T^ lying in the range of 4-10 MeV are listed under the labels Si and ^2 respectively. 
For the case of ^^^Sm the transition to the state at 4.08 MeV is included in Si (see text). 

The distribution of the orbital strength predicted for the QRPA states lying below the limit of 
4 MeV is compared with the experimental data in Fig. 6. Several features revealed by this picture 
should be mentioned. The number of states shown by experiment is equal to the corresponding 
theoretical number. Exception is ^^^Sm where the first two small experimental sticks correspond 
to one relatively big stick. Even in this case, the dominant 4 peaks bunched around 3 MeV are 
well simulated by the predicted peaks in the range of 2.9-3.6 MeV. Although the main peaks are 
shifted a little with respect to the corresponding experimental results, the shape of the theoretical 
and experimental strength distribution are in good agreement with each other. Also we notice that 
going from the deformed nucleus ^^^Sm toward the spherical case of ^^^Sm, the orbital strength 
is less fragmented which is consistent with the fact that deformation causes fragmentation. The 
transition probabilities to all states from the interval 2-10 MeV are shown in Fig. 7. While in -"^^^Sm 
the strength is almost continuously distributed and shaped in maxima alternating with minima, 
in the near spherical nuclei the interval containing strength is decreased and moreover the groups 
of peaks are well separated and less members are included in a group. Several doublet structures 
are noticed in the lightest two isotopes. From our analysis it results that one cannot state that 
there is a sharp separation between orbital and spin like excitation and that the border is 4 MeV. 
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On the contrary, one could say that for deformed case there are reasonable chances to find orbital 
transition beyond 4 MeV while for lighter isotopes the spin transition could be visible below 4 
MeV. For example in ^^^Sm the transitions to the states of energy equal to 3.69 and 4.03 MeV 
are of proton spin nature while the state at 5.419 MeV is of orbital nature. Indeed for the latest 
transition the ratio orbital to spin is 2.11 and its magnitude is 0.33 hn^. In ^^^Sm the states 
with the energies 2.79 and 3.64 MeV have spin-neutron and spin-proton character. ^^°Sm has 
also spin-flip like states below 4 MeV, at 3.05 MeV, and an orbital state of energy equal to 4.005 
MeV. In our formalism none of the states below 4 MeV calculated for ^^^'^^^Sm is of pure orbital 
nature. Indeed, for ^^^Sm the two states lying below 4 Mev in Fig 7 are characterized by the ratio 
Ro/s = 0.54; 0.26 while for ^^^Sm this ratio is equal to 0.324. In all three cases the spin contribution 
is brought by protons. The highest peak in the heaviest three isotopes is of neutron nature and is 
lying in the interval 7.5-8MeV. The second highest peak in ^^^Sm is dominated by the proton spin 
component but has also neutron admixture. Its energy is smaller than that of the neutron peak, 
by 700keV. In the next two isotopes there is no second peak of height comparable to that of the 
neutron peak. In ^^^Sm there are two proton peaks with energies 4.05 and 7.51 MeV with strength 
slightly larger than 0.5 fiN^ and a neutron peak at 5.5 MeV with a strength of about 0.5 /j-n'^. The 
lightest isotope exhibits a dominant neutron peak at 6.35 MeV, two proton peaks lying close to 4.5 
and 9 MeV and three neutron peaks with energies 5.3, 5.4 and 6.6 MeV. Our calculations suggest 
that one cannot separate sharply an energy interval where only neutron spin strength can be found 
and one only with proton transitions. However a certain polarization of the strength distribution 
may be seen. For example, for the first three nuclei (^^^'^^'^^^^^Sm) beyond 7.5 MeV transitions are, 
with very few exceptions of small proton transitions, of neutron type. In the last two ( ^^^'^^^Sm) 
cases in the interval mentioned above, the situation is opposite. An interesting feature one sees in 
^^^Sm is that the regions with proton and neutron like transitions come in successive order. Indeed 
in the interval 4-5 MeV one observes only one proton peak. The next interval 5-7MeV is populated 
mainly by neutron transitions, the protons appear to be impurities being the smaller components 
in the last two doublets. Finally, beyond 7 MeV there are only proton like transitions. 

One of the most interesting properties of the scissors like excitation for Sm isotopes was pointed 
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out in Ref. Q, and refers to the linear dependence of the total strength below 4 MeV, on the 
deformation squared. This dependence is nicely obtained in Fig. 8 and this suggests once more a 
very good agreement with the experiment. Such a relationship between the total orbital strength 
and nuclear deformation is a rigorous prove that the magnetic states populated in the inelastic 
(e,e') backward scattering are of collective nature. On the other hand, the E2 transition probability 
from the ground state to the first collective 2"*" state has also a linear dependence on the nuclear 
deformation squared. This suggests a proportionality relationship between the two quantities. 
Indeed in Ref. (5^j5^ it has been shown that the proportionality factor is the average energy 
defined as 

=Y,EiB,{Ml,^)/J2Bi{Ml,^). (4.3) 

The theoretical and experimental average energies, defined as above, are compared with each other 
in Fig. 9. The agreement is again very good. Indeed the deviation of theory from experimental 
data vary from lOOkeV in the case of ^^^Sm, to 210 keV in ^^"^Sm. The collective features of the 
scissors states might be also evidenced in the plot of the total orbital strength below 4 MeV versus 
the factor 



introduced by Casten in Ref. Indeed, as shown in Ref. [^J] the graph obtained is almost 

identical to that representing the BE2 values as function of the same quantity. In the above 
equation A'p, A^^ denote the proton and neutron boson numbers respectively, defined in the IB A 
model as being equal to half the nucleon number in the open shell. In our case the predicted and 
experimental values were interpolated as in the reference quoted above by a continuous curve. We 
note that here the saturation process start with ^^^Sm at P ~ 5 and the limit strength is 2.9 /u^. 



V. CONCLUSIONS 



A model Hamiltonian has been treated within the QRPA formalism based on a projected 
spherical single particle basis. The mean field for the single particle motion is defined by averaging 
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a particle-core Hamiltonian over a coherent state describing a deformed phenomenological core. 
The monopole part of the couphng term is determined from the volume conservation condition. 

An orthogonal basis is defined through projection technique from a product function of a spher- 
ical single particle state and a coherent state for the quadrupole shape coordinate describing the 
phenomenological core. Although this basis comprises both the single particle and core coordi- 
nates, it may be used for treating the operators acting only on particle coordinates by integrating 
on the core coordinates. A set of single particle energies are defined by averaging the mean field 
Hamiltonian on each member of the basis. These energies depend on deformation in a similar way 
as the Nilsson energies. It is argued that these energies approximate quite well the eigenvalues of 
the mean field Hamiltonian while the projected states are close to the states which are obtained 
by projecting the good angular momentum from the eigenstates of the mean field Hamiltonian. 

The QRPA formalism provides the dipole magnetic states. Special attention is focused to those 
states lying in the interval 0-10 MeV. For these states the B(M1) values were calculated by using 
the transition operator given by Eq. (3.5). 

Applications were made for the even-even isotopes of Sm. The intervals 2-4 and 4-10 MeV 
were considered separately since the states in these intervals are excited by different parts of the 
transition operator. 

For the first interval we compared the predictions with the corresponding data for the following 
observables. The distribution of the theoretical orbital strength has similar shape as the experi- 
mental results. Moreover the number of the transitions seen experimentally is equal to the number 
of transitions of a visible size in the plots shown. 

The very good quality of the agreement between theory and experiment can be seen from the 
plot showing the quadratic dependence on the nuclear deformation of the total orbital strength. 
The same conclusion can be drawn from the plots representing the average energy in the interval 
0-4 MeV as well as from the graph of total strength as function of the Casten's factor P. 

Concerning the spin transitions we considered first the case of ^^^Sm where experimental results 
for strength distributions are available. As shown in Fig. 5, the predicted distribution agrees quite 
well with the experimental one. For all considered isotopes the spin transitions are presented in 
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Fig. 7. One clearly notices how the fragmentation is influenced by deformation and moreover 
the appearance of the shell structure when one goes from the deformed to the spherical isotopes. 
Comments about the isotopic dependence of the spin strength and upon the fact whether there are 
orbital transitions beyond 4 MeV or spin transitions below 4 MeV are included. 

We stress that the formalism used in the present paper is similar to that from Ref. except 
for the single particle energies. Indeed there they have a linear dependence on deformation while 
here higher powers of deformation are included by implementing the volume conservation condition 
for the mean field. This feature seems to be essential in obtaining a better description of both the 
orbital and spin strengths. 

As a final conclusion one may say that the present model, using a projected spherical single 
particle basis and a mean field defining a set of deformed single particle energies which depend 
non-linearly on deformation, is an efficient method for a realistic description of the Ml properties 
caused by convection currents and spin distributions in nuclei. 
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FIG. 1. Proton single particle energies given in units of Hloq. 
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FIG. 2. Neutron single particle energies given in units of Hloq. 
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FIG. 3. Neutron quasiparticle energies given in units of MeV. 
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FIG. 4. Neutron quasiparticle energies given in units of MeV. 
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FIG. 5. The theoretical histogram of the spin strength (lower panel) is compared with 
experimental strength distribution. The experimental plot is taken from Ref.[9] 
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FIG. 6. The predicted orbital strength (lower part) is compared with the experimental data 
(upper part). Experimental data are from Refs[2-9] 
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FIG. 7. Theoretical Ml strengths, in the range of 2 to 10 MeV, are plotted as function of energy 
for some even-even Sm isotopes. 
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FIG. 8. The orbital total strength in the range 0-4 MeV is presented as function of 5^, for the 
even-even Sm isotopes. Both the theoretical and experimental results (Refs. [2-8]) are shown. 
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FIG. 9. Theoretical (triangles) and experimental (black circles) average energies calculated by 
means of Eq. (4.3)are given for the 5 isotopes of Sm, studied in the present paper. 
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FIG. 10. The orbital Ml strength in the interval 0-4 MeV, as function of the ratio P defined 
by Eq. (4.4). 
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